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We calculate the non-forward quark matrix elements for operators with two covariant derivatives in one-loop 
lattice perturbation theory using Wilson fermions. These matrix elements are needed in the renormalisation 
of the second moment of generalised parton distributions measured in lattice QCD. For some commonly used 
representations of the hypercubic group we determine the sets of all mixing operators and find the matrices of 
mixing and renormalisation factors. 



In recent years generalised parton distributions 
(GPDs) have been intensively studied both in 
experiments and theoretically pp to extend our 
knowledge of hadron structure. However, the di- 
rect experimental access to GPDs beyond the lim- 
iting cases of distribution functions and simple 
form factors is limited. Therefore, it is indispens- 
able to obtain complementary information from 
lattice QCD calculations. 

On the lattice we can compute matrix ele- 
ments of local composite operators, and mo- 
ments of GPDs can be related to such matrix el- 
ements taken between states of different nucleon 
momenta and spins. First results for moments 
of GPDs have been published recently see 
also 0, and soon results from improved calcu- 
lations should become available. 

In order to relate lattice measurements to con- 
tinuum quantities we have to investigate the 
renormalisation and mixing of the operators in- 
volved. When non-forward matrix elements are 
studied, new features arise, which make a recon- 
sideration of the renormalisation problem neces- 
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sary. In particular, the mixing with "external 
ordinary derivatives" , i.e. with operators of the 
form d^dv • • ■ (ip ■ • ■ ip), needs to be taken into ac- 
count. On the lattice the mixing patterns are 
usually more complicated than in the continuum, 
because covariance under the hypercubic group 
H(4) imposes less stringent restrictions than 0(4) 
covariance. The necessity to consider also exter- 
nal ordinary derivatives enlarges the set of con- 
tributing operators even further. These compli- 
cations do not yet arise for operators with zero 
or one covariant derivative. Hence for these mo- 
ments the renormalisation factors can be taken 
over from the forward case. 

Here we investigate the renormalisation prob- 
lem for operators with two covariant derivatives 
sandwiched between off-shell quark states at dif- 
ferent momenta within the framework of one- loop 
lattice perturbation theory and report on first re- 
sults 0]. We use the Wilson plaquette gauge ac- 
tion with Wilson fermions and perform the cal- 
culations in Feynman gauge. 

Let us first discuss renormalisation and mixing 
in general. Denote by Tj(p',p, fi,gn,e) (with j = 
1,2,..., N) the dimensionally regularised ampu- 
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tated vertex functions of N mixing operators Oj , 
p and p' are the quark momenta. The correspond- 
ing Born terms are denoted by Tj(p',p). The 
renormalised coupling constant g^ is related to 
the bare coupling g by g^ = fi~ 2e g 2 (l + 0(g 2 )), 
where is the renormalisation scale. In the 
MS scheme the corresponding renormalised ver- 
tex functions are given by 

N 

Tf(p', P ^,g R ) = rf (p', P ) + 9 i\ Ysi-'yjk) 

L " (1) 



k=l 



x\n ( ^Lr*(p>,p) + f 3 (pi,p)], 

where fj is the finite one-loop contribution. All 
O(g^) contributions are neglected. 

In the absence of mixing with lower- 
dimensional operators the vertex functions reg- 
ularised on a lattice with lattice spacing a are 

N 

Tf(p',p,a,g R ) =r»)+. 9 |[5> 7 ^) 

fe=i ( 2 ) 
xln(ia 2 (p'+p) 2 )rf(p',p) + //(p',p) . 

The relation between T k and the MS renor- 
malised vertex functions Tf is given by 
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where the matrix £ is found to be 

Cjk = lj k In (a 2 (i 2 ) - cj k 

with the constants cj k determined from 
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(3) 
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(5) 



f J L (p',p)-f j (p',p) = ^cY k T k i (p',p). 
fe=i 

Mixing with lower-dimensional operators leads to 
the appearance of additional terms on the r.h.s. 
of Eq. ©, for details see gj. 

The operators potentially contributing to the 
mixing have to transform identically according to 
a given irreducible representation of 0(4) or H (4) 
and should have the same charge conjugation par- 
ity. 

In terms of the quark wave function renormali- 
sation constant Z^ and the matrix Zj k of mixing 
and renormalisation coefficients, the connection 
between the bare lattice vertex functions and the 



MS renormalised vertex functions can be written 
as 
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Using the known Z., 



Z% = \ — 



9lC, 



16tt 



f [In (a V) + 1 + ox] 



(6) 
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with <7l = 11.8524 for Wilson fermions, we find 
9rC f 



Zjk 



16tt 2 



[jjk log(a 2 ^ 2 ) + cjk] 



(8) 



with jjk = Sj k - (167r 2 / 'C 'F)jJ k and c jk = 5 jk (l + 
a L ) + (167r 2 /C F )cJ k . 

The one-loop computation has been performed 
symbolically, adopting and significantly extend- 
ing a Mathematica program package developed 
originally for the case of moments of structure 
functions using Wilson p], clover |S] and over- 
lap fermions 0. Using that approach we have 
full analytic control over pole cancellation. The 
Lorentz index structure of the matrix elements is 
left completely free, so that we are able to con- 
struct all representations of H(A) for the second 
moments in the non- for ward case. 

We have to deal with operators of the following 

general types (D=D - D)- 



0°Z = (-V*)j>^D v D„iP, 

o?,i, = 



O 



flV 'U) 



= (-1/4)9^(^7^) , 
= (-i/2)flu,(#y/*,7„ty) 



(9) 



and similarly for operators with 7^75. In lattice 
momentum space we realise the operators with 
non-zero momentum transfer q (showing as an ex- 
ample an operator with one covariant derivative) 
through 

(i>D^){ q )= l -Y. ciq{x+aN2) * (10) 

[ip(x)U XtlJ/ ip(x + ajj) - tfj(x + afi)Ul^ip(x)] . 

Eq. I|10l) basically defines the Feynman rules in 
lattice perturbation theory. 

Consider now the following operator: 



Cl — Cfn4} _ 7J ( C{224} 
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Cjfc 



V 



-12.1274 


3.3060 






1.4913 
20.6178 
-8.01456 


3.2644 
3.2644 



Table 1 

The finite mixing part for the operators Oi, 



0.3685 


-14.8516 






-0.4160 


4.3023 
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0U, B3J- 



0.0156 
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-0.9285 


0.3501 
0.0050 



0.1498 \ 


0.7380 


0.0149 
0.3600 / 



with charge conjugation parity C = — 1 and cor- 



responding to the representation r 
use the notations 151 : 
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of H{4). We 
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The following operators (transforming identi- 
cally) mix in one- loop order with 0\ : 



2 - C{fi4} 
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Or = O 



y 5 — ^1121311 ' ^ 6 ~~ ^<<213)> 

together with the lower-dimensional operator 



O s = O c A 



o: 



I) 



O a 433 ) /2. 



(14) 



'411 V^422 ~ ^433; 

An additional operator with the same symmetry 
properties does not contribute in the one-loop ap- 
proximation. The anomalous dimension matrix is 
obtained as 
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and the finite mixing matrix part is given in Ta- 
ble ^ There is an additional mixing of 0\ with 
the lower-dimensional operator Os (in Born ap- 
proximation) : 



01 



-part 



1^,-0.5177, i | 



(16) 



This mixing leads to a contribution which poten- 
tially diverges like 1/a in the continuum limit. 
The operator Os has to be subtracted non- 
perturbatively from the matrix element of the op- 
erator 0\ which might be a difficult task in sim- 
ulations. 

In summary, we have presented first results for 

the one-loop quark matrix elements of operators 

needed for the second moments of GPDs. This 

allows us to determine the mixing matrix of Z- 

factors in the MS-scheme for operators commonly 

used in simulations. In the case of the represen- 
ts 1 ) 

tation t\ , a sizeable mixing with operators of 
the same dimension appears. Moreover, mixing 
with a lower-dimensional operator has to be taken 
into account. Additional results, in particular for 
transversity operators, will be presented in 0]. 

This work has been supported by DFG un- 
der contract FOR 465 (Forschergruppe Gitter- 
Hadronen-Phanomenologie) . 
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